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Energy transport mechanisms can be generated by imposing relations between null 
tetrad Ricci components. Several kinds of mass and density transport generated by 
these relations are studied for the generalized Vaidya system. 


INTRODUCTION 


Q 


The Vaidya |l[ spacetime generalized the Schwarzschild vacuum solution by allowing mass 
parameter mo to be a function of retarded time. The extension created a spherically symmet¬ 


ric null fluid atmosphere. Glass and Krisch [2j pointed out that allowing the mass function 
to also depend on the radial coordinate created a more complex atmosphere containing an 
anisotropic string fluid in addition to Vaidya radiation. 

The fluid parameters for the Vaidya atmosphere {p, p r , p±} all depend on time through 
their relationship to the time dependent mass function m{u,r ); they control the motion 
of matter through the atmosphere. The Ricci tensor null tetrad components (i.e. Ricci 
spinor components), <f>n and $ 22 , for the generalized Vaidya spacetime also depend on 
both spatial and time derivatives of the mass function. If we impose relations between these 
components then those relations will, in turn, impose matter transport in the Vaidya system. 
Conversely, if a particular mode of matter transport were assumed for the mass or one of 
the fluid parameters, it would contain within it a Ricci relation. Relating Ricci components 
is a way of imposing and classifying a range of atmospheric matter transport mechanisms. 
In this paper we consider the generalized Vaidya metric ^jsj], and examine how functional 
relations between $ n and $22 lead to mass transport described by the diffusion equation, 
by the wave equation, and by a dissipative transport equation of Telegrapher type. 

In studying fluid transport, the transport equations are often reducible to ordinary dif¬ 
ferential equations by introducing a similarity variable; a well known example of this is 


the diffusion variable, ~ r/\/t , suggested by Boltzmann 0] in 1894. All the transport 
equations studied in this paper have similarity solutions. The physical similarity of matter 
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. The matter transport considered 
here is in the atmosphere around a compact object or an already existing black hole. By 
examining the similarity structure of the mass solution to the transport equations, we are 
able to relate the function of proportionality between the Ricci components to the spatial 
part of a similarity variable. 

The paper is structured as follows; In the next section we briefly review the generalized 
Vaidya spacetime. Matter transport mechanisms and their similarity structure are discussed 
in section III. Some explicit examples are given in section IV. Metric and tetrad details are 
provided in Appendix A 

Our sign conventions are 2 A v . t y a p\ = A^R^ ua/3 , = R a ^ va , and metric signature 

Greek indices range over (0,1,2,3) = (u, r, d, </?). m abbreviates dm/du, m' 
abbreviates dm/dr, with another prime for each higher derivative. Overhead carets denote 
unit vectors. Field equations are G /JiV = SnT^. 


fields has been discussed by many authors 


5 ], 6 ], 7 ], 


II. GENERALIZED VAIDYA SPACETIME 


The Vaidya metric generalizes vacuum Schwarzschild when m 0 —)► m{u) 
g v a fdx a dxP = A vad du 2 + 2 dudr - r 2 (dd 2 + sin 2 Hdip 2 ) 


( 1 ) 


where A vad = 1—2 m{u)/r. The only non-zero Weyl null tetrad component is \&2 = —m(u)/r 3 . 
The Ricci tensor, with l a dx a = du, is 


pvad _ 7 7 

-0.^/3 o fcf yl> 

r 


2m , 


■a/3 ~~ ^2 ''“O' 


The Vaidya metric has been generalized so that m{u) —> m{u,r ) and A vad —> 
A^ = 1 — 2 m(u,r)/r. 

g///pdx a dx 0 = A^du 2 + 2 dudr — r 2 (dd 2 + sin 2 dd(p 2 ). (2) 

The metric is Petrov type D with l M and n M principal null vectors (details are in Appendix A). 


l^dx^ 1 = du, 

(3a) 

n^dx 11 = (A gv / 2) du + dr, 

(3b) 

m^dx^ = —(r/y/2)(dd + i sin ddip). 

(3c) 
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Basis vectors for the generalized Vaidya metric are the unit vector set (v^, r M , tp^) and 
the related null set (/ M , n M , m M , m M ) such that 

gtjZ = vk - f^ v - ’&pd v - (p ll (p v (4) 

= Ifjiu + n^lv - m^m v - rh^m u . 


In terms of the basis vectors, metric has energy-momentum 



-87t Tj% = W*' + + p r r^r u + + (p^v) 

(5) 

with components 

47T0 = —777,/r 2 , 

(6a) 


4^ 

-1 

1 

■i 

Ah 

LO 

(6b) 


8np± = —m"/r. 

(6c) 

The Einstein tensor is computed from g^ and given by 



G% = -2$ n (^n„ + n^l v + m^m v + m^m v ) 

(7) 


- 2$ 22 l/il v ~ (Jl/1 )g%, 


with components 

$11 = (2 m! — rm")/(4r 2 ) = 2ir(pj_ — p r ), 

(8a) 


$22 = -m/r 2 , 

(8b) 


IZ = 2(rm" + 2m!) /r 2 = —167r(p_L +p r )- 

(8c) 


An inspection of Eq.(JHJ) shows that a relation between $n and $22 will generate matter 
transport. 

Sfj, u = ~ \'^' SV 9^u) is the trace-free Ricci tensor. Its eigenspectrum determines the 

Segre type of metric gf ^. The characteristic equation, det[S' — XI] = 0, is 

A 4 — ( < f ) 22)A 3 — 2(<h 2 1 )A 2 + ($? 1 < h 2 2)A + = 0, (9) 

with factors 

(A — $n)(A + <f>ii)(A 2 — <f >22 A — $ 1 !). (10) 

This set corresponds to Segre type [2,(11)]. The Ricci relations established below will not 
change the Segre type. 
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III. GEOMETRY AND TRANSPORT OF MATTER 


A. Mass diffusion from fluid relations 


The fluid components of the energy-momentum in equation (JHJ) depend on both time and 
spatial derivatives, so matter transport will occur. Implicit in these equations is the relation 

47rr 2 p = d r m. 


If one also assumes 


rh = 4:iiDor 2 d r p 


then both the mass and density obey diffusion equations, the mass diffusing in a space whose 
determinant is dual to Vaidya jij. 


p = D 0 V 2 p, (11a) 

m = D 0 'V~ 2 m (Hb) 


with V 2 = r 2 (d / dr)r 2 (d / dr ), V 2 = r 2 (d/dr)r 2 (d/dr), and D 0 the diffusion constant. 
Equating the non-zero Ricci components will also produce motion of the mass and density. 


B. Mass diffusion from Ricci components 


Consider the quantities (dr^n) and (r 2 *!^)- Equations (IKal) and 
tions) allow the relation between these quantities to be written as 

<9 r (r 2 $- 


0u(4r$n) = r~d r [- 


■22 J 


(with no assump- 


( 12 ) 


If there is a general linear relation between and $22 such as 

^22 = hd(r ) 4 >n. 


(13) 


then both Ricci components will evolve diffusively. We will show in the next section that 
hd(r) is related to a spatial similarity variable for the diffusion equation. For exact solutions 
mass is a more physical quantity and, if we impose this relation and examine mass transport, 
equations (ISafl and (l8b|) imply 

<f>n = —{r/4)d r {m/r 2 ) 

<£>22 = — m/r 2 = — {rhd/4)d r (rn!/r 2 ). 


(14a) 

(14b) 
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The linear relation m together with (I14M yields 

rh = (r 3 hd/4:)d r (m' / r 2 ) . (15) 

and so the mass will also move diffusively. 

The homogeneous solution of Eq. csD is 

m = r 3 ci(u) + c 2 (w), c 2 + r 3 ci = 0 (16) 

8ttp = —8-Kpr = 6ci(w) 
pr = Pa- 

Equation c 2 + r 3 ci = 0 is satished only for c\ and c 2 both constant. 


C. Mass transport by wave motion 


The basic relation between the Ricci components, Eq. (ED, is 


d n (4r<f>ii) 


r 2 d r 


.<9 r (r 2< f) 


22 J 


If one assumes the Ricci components are related by 


d u &2-2 — h w (r )&ii (17) 

then r 2< h 22 will evolve as a wave equation solution. Focusing again on the mass behavior, 
Eq.(JHH) for <hn and Eq. (J8bl) for $ 22 provide a wave equation for the mass: 

bl — ( h w /4)r 3 S7~ 2 m = 0 


or 

rh — (h w /A)r 5 d r (r~ 2 d r m) = 0. (18) 

Classical wave motion provides shape preserving traveling solutions to the wave equations 
for functions with argument (kr — u:t). Therefore, consider a variable of the form r] = 
R(r ) —T(u). The wave equation becomes 

[T 2 - (A„/4 )r 3 (R') 2 ]m m — [T + (h u ,/4)(r 3 R" - 2r 2 R')]m„ = 0. 


(19) 



6 


The transport equation is 


oim wl — /3m v = 0 
T 2 - (h w /A)r\R') 2 = a 
f + (h w /4)(r 3 R" - 2r 2 R!) = (3. 

The simplest solutions require a and (3 to be separately zero. These are 

h w = r 3 ^) 2 ’ = C 2' r3 + c 3 ) T = T 0 u + T 1: (20) 

and, as in diffusive transport, the proportionality function between the Ricci components 
is related to the similarity variable. The mass is given by any function of argument 77 , 
m = F(r)). 


D. Dissipative transport 


It is clear that the mode of matter transport depends on the relation chosen between 
$22 and $ 11 . Dissipative waves obey a Telegrapher equation P|, combining both wave and 
diffusive elements. 

X - CqV 2 x + cix = 0 (21) 

where V 2 x = (l/r 2 )d r (r 2 d r x)- The relationship 

{d u + ao)<h 22 + ht(r)$ 11 = 0 ( 22 ) 

will generate an equation of Telegrapher form for mass 

m — {h t /A)rS/~ 2 m + a 0 m — 0. (23) 


This form of causal dissipative transport provides a richer tool for modeling than a pure 
diffusion equation, since it avoids the infinite propagation speeds associated with parabolic 


diffusion equations 


0 . 


0 - 


Causal transport has recently been discussed by Herrera and Santos 


As before, the function of proportionality h t is related to a similarity variable. 




7 


E. Similarity and hd 


To see the relation of the proportionality function to the similarity structure of the dif¬ 
fusion equation, assume a mass solution and similarity variable with forms 

m — rn 0 = T(u)F(r /) 
i] = R(r)/L(u). 

Rewriting the mass diffusion equation (USD, with F v = dF/dri results in 


vv 


R" 2 1 AR 
~ Tr* + h d r(R') 2 


This diffusion equation should have similarity form 


0 


>T 


hd'r(R') 2 T 


F = 0. 


(24) 


F m + 2arjF v - j3F = 0. 


We look for solutions with the same similarity form, and so the coefficients impose con¬ 
straints. The coefficient of F relates (3 and hd(r) 


q — _r 2_ 

; h d r(R') 2 T' 

Removing r-dependence provides the relation between R[r) and hd{r ) : 

/ - 4c o 

d ~ r(R ') 2 ' 

With this hd constraint, the u-dependent part of the coefficient (3 implies 


(25) 


(26) 



(27) 


This allows Eq. (El to be written as 

R" 


F 


vv 


L{ 


(R'Y 


2 

rR' 


} + (R/co)L 


— (3F 


0 . 


The coefficient of F v requires 

2a,, = 2 a(R/L) = L[-F- - -|_] + (R/cJL. (28) 

Solutions of this equation depend on the value of a. Some example solutions will be given 
in Section IV. We note here that the choice R = R^r 3 + R \, which solves 


R"/R’ - 2/r = 0 , 


( 29 ) 















provides an analog of the Boltzman similarity variable. For this choice Eq. (I28|) yields 


with similarity variable R(r)/L{u) 


and proportionality function 


— Acoau T Lt* 


V = 


R 0 r 3 + R,\ 
\J 4c 0 cm + Ll 

4 c 0 


hd = 


r 5 Rl' 


If Co = D 0 Rl then hd = 4 D 0 /r 5 . This choice yields a simple diffusion equation for mass. 


F. Similarity and ht 

Assuming a mass and similarity variable of the form 

m — m 0 = T(u)F(r]), rj = R{r)/L{u) 


we find 




T v „ ,h t . 


,ht 


/\2 


)r(R') 




R(2L— +L- 2— + a 0 L) + (-f )L(rR" - 2 R’) 

1 Lj 4 


T 


T, 


L (— + «q — ) 


F = 0 


'T T' 

The relationship of h t to the scaling variable in this similarity equation depends on the form 
of the time scaling function. The example that we shall use in the next section has L = L 0 , 
and for this choice the similarity equation becomes 


-jrm 2 F„ + - rR")F„ + ig(T + ao l)F = 0, 


(30) 


For this choice of L, we see the diffusive relation between the proportionality function and 
the spatial scaling function again emerges as 

h ‘ = < 31) 

Telegrapher transport could allow many other relations to be imposed, reflecting the 
richer solution structure of this transport mechanism. In the next section we give some 
examples. 
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IV. TRANSPORT EXAMPLES 

A. Similarity Solutions 

As an example of the mass transport solutions based on similarity, we consider a single 
scenario to which we apply all three transport mechanisms. The physical setting considered 
is a compact object of mass m(u, rj) with an atmosphere in which transport is occurring. We 
use solutions where the scale variable R and proportionality function h(r) are the same for 
all transports. 

fl=A- s , L — L 0 ft(r) = AL (32) 

3 /to r 

In simple diffusion problems, the choice L = Lq is applied to bounded systems where Lq can 
be identified with a natural physical scale. For the Vaidya black hole, the Schwarzschild 
radius provides a physical distance scale and we could identify L 0 with an associated distance. 

1. Diffusive transport 

For diffusive transport with the choices above, the equations describing the evolution of 
the mass function are 

F m — (3F — 0, T/T = (3c 0 /L 2 0 . (33) 

The mass function, with time parameter rf 1 = /3D 0 Rq/ Lq, will be given by 

m = m 0 + T 0 e U ^ T ° F^rf) (34) 

F{v) = F o sm( \fW\ r l + 5), P < 0 

= F 0 sinh(A /fir) + 5), /3 > 0. 

For /3 < 0 , the atmosphere is decaying in time as the Vaidya photons carry energy out of 
the system, while for /3 > 0, mass is accreting. From the held equations, the density of the 
atmosphere is described by 


p = Poe U,T ° cos(a/[5[?7 + 5) 
= Poe u,TO cos1i(a /Pp + 5) 


( 35 ) 
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where we have identified the initial atmospheric density and a time constant for the accretion 
or decay 


Po = 


T 0 F 0 ^/W\ 


2n Ln 


(36) 


The similarity variable in this example is simply a distance coordinate. The decaying 
solution could describe a bounded atmosphere whose density decreases with distance away 
from the surface. There is zero radial pressure outer boundary described by \J\fi\\ t]b + £ = 
7 t/2 . In the decaying solutions, the atmosphere will go asymptotically to zero leaving a 
compact object of mass mo- Because of the outgoing Vaidya radiation, there is no vacuum 
match until the atmosphere is gone. For the accreting solution, there is no zero radial 
pressure boundary. Since the accreting mass is entering the atmosphere from the exterior, 
this is expected. The density profile depends strongly on the phase 5. For <5 = 0, the density 
increases as one looks upward from the core surface. For non-zero 5, the density decreases 
going up from the surface, reaching a minimum value which might be identified with the 
boundary between the atmosphere and the source of the accreting mass. 


2. Telegrapher transport 


The equations for Telegrapher transport with the scaling functions chosen are 


L 


2 

0 


Co 


F m - (3F = 0 

T T 
(— + a 0 —) = (5 


Telegrapher transport has a wider set of solutions for the scaling choice than pure diffusive 
transport. One solution, using (i — 0, is 


T = T 0 e~ aou (37) 

F = F 0 rj 

m = mo + pi-r 3 e _ “ 0U , pi = T 0 F 0 /L 0 . 

The atmospheric density is related to the mass through the field equations by p = m! /( 47 rr 2 ). 
The density for this case does not vary with radial distance from the black hole but does 
reflect the atmospheric decay, approaching zero as the atmosphere vanishes and the mass 
becomes m 0 : 


p = p 0 e Po = 3-Rop,/4jr. 


(38) 
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The /3 7 ^ 0 solutions are similar to the diffusive solutions 

T = T 0 e 7 “, P=(L 2 0 / Cl ){7 2 + «o7 ) (39) 

F = F 0 sin {\fW\n + 5), P < 0 
F = F 0 sinh( \fjdr) + 5), P > 0 
m = m 0 + T(u)F(rj). 

B. A diffusively evaporating atmosphere 

a. The boundary behavior ♦ In this example, we consider a compact object whose 
atmosphere is diffusively evaporating. An exact mass solution, given in [2], is 

m = m 0 + (4n/3)r 3 p 0 — (4:71 / 3) k 2 (r 5 /10 + D 0 ur 3 ). (40) 

For this example the proportionality function follows directly from the diffusion equation 

h d = 4D 0 /r. 

The atmospheric density associated with this mass is 

p = po - k 2 (r 2 /6+ Q 0 u). (41) 

A boundary can be defined by requiring the radial pressure (and density, from the equa¬ 
tion of state) to be zero. Solving the density equation for a boundary radius gives the 
boundary as a function of retarded time 

6 

Rl = — (po - k 2 D 0 u). 
k 2 

Let the boundary move inward as the atmosphere evaporates until the atmosphere is gone 
and the radius is at the compact object boundary, R c . This happens in time uq, thus 

Rl = -j~(po ~ k 2 T>QUo). (42) 

k 2 

It follows that 

Rl(u) = R 2 c + 6 D 0 (uo - u). (43) 

The bounding surface of the core is parameterized as R c = 2mo«. Substituting into (H3J) 
provides 


Rl(u) = 4mla 2 + 62) 0 (u 0 — u ). 


( 44 ) 
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The mass function places a constraint on the parameters. From Ea. (1421) at time Uq 


/ u 128, 5 5 

(a — 1 )m 0 = -^rk 2 nm 0 a . 

Using Eas-fl-Hl and (TT31) . the time for the evaporation process to complete i 


D 0 m 0 = 


647rp 0 mQa 3 

_ 15(a- 1) 


— 1 


The density and mass evolve as 


ko 


P = t r(R c ~ R b ) + k 2 T>o(uo - u ) 
o 


(45) 

(46) 

(47) 

45 "' 0 ' 3 c ^ 

b. Extremal case ♦ The time at which the evaporation of the atmosphere is complete 

has an extremal value. Extremising Eq. (ES1) one Ends 

2 (a — l) 2 647rp 0 ?7io 

a 3 (4a — 5) 15 

The extremal time for the evaporation is 

2 


m = mo + 




' rp2 _ d2 

^ + Do(n 0 -«) 

6 


The second derivative shows this is a minimum. For positive times, we require 


(49) 


1.25 < a < 1.5 


c. Time estimates ♦ Take m 0 to be a solar size : mo ~ 1.5 x 10 3 m. The diffusivity 
constant 33o is related to the constancy of the jump distance (L) and frequency (/) 4], D 0 
can be estimated as 


Sr 


L 2 f. 


For example, if the Debye frequency in a solid is the same as the jump frequency, we would 
have / ~ 10 13 Hz and the jump distance could be of the order of an Angstrom, so that an 
estimate for the diffusivity is 

33o ~ 10~ 7 m 2 /s. 

Using these estimates, the time for the atmosphere to vanish by diffusion is roughly 

Mo = 10 7 (ba 2 x 2.25 x 10 8 [f^] 

3 4a — 5 

1.5 x 10 13 [^^]« 2 

4ck — 5 


rs./ 
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For a = 1.4, the time is about 10 13 sec ~ 10 5 years. The diffusivity is much larger than is 
normally measured since not all atomic oscillations will have an associated jump. A smaller 
diffusivity would increase the time. 


V. DISCUSSION 


In this paper we have examined some of the atmospheric matter transport mechanisms 
introduced by imposing a relation between Ricci components of the generalized Vaidya 
spacetime. We found that diffusive, wave, and general Telegrapher-type transport can all 
follow from such a relation. Imposing a Ricci relation is a way of unifying matter transport 
mechanisms; if one started with a particular transport mechanism then a Ricci relation 
would emerge. 

Examples of the matter transport mechanisms were given. All of the exact solutions 
describe atmospheres that are either accreting or decaying. The first two examples had 
the same similarity function but the density profiles in each solution were very different. 
The first Telegrapher solution describes a core object with a decaying atmospheric density 
varying only with time, and whose tangential and radial stress are all spatially constant 
tensions. In the pure diffusive solution, the densities vary with both distance from the core 
surface and time. The decaying atmospheres could be bounded by an exterior fluid of Vaidya 
radiation. The accreting solutions were surrounded by the source of accreting matter and 
Vaidya radiation. The third example, with a different proportionality function, examined a 
diffusing atmosphere around a core object of radius R c in the range 2.5m 0 < R c < 3 m 0 . 
For an R c = 2.8 m 0 object with a diffusivity based on a solid Debye frequency there was 
a minimum evaporation time of about 10 5 years. Jump frequencies based on the actual 
atmospheric fluid would increase this time. Apart from the actual size of the evaporation 
time for a specific object, the model predicts that atmospheres around smaller core objects, 
near the lower end of the range, will take very much longer to diffusively evaporate than the 
atmospheres of larger core objects. This can be understood by looking at the initial size of 
the atmosphere. Using Ea. dlUl at u — 0, Eq. flUJ), and the parametrized boundary surface 
R c = 2mo«, the initial atmospheric radius is 


= 2 R 


2 

c 


a — 1 
4a — 5 


Rb |u=0 


1.25 < a < 1.5. 
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The smaller core objects have larger atmospheric envelopes, taking longer to evaporate. The 
smaller core has an extended atmosphere because its gravitational field is not strong enough 
to hold the atmosphere compactly. 

Some additional insights into the meaning of the proportionality function can be seen for 
diffusive and wave transport by writing the transport equations in terms of variable w = r 3 . 
The diffusive and wave transport equations become 

,9 r 5 h d 

m = )m ww 

Qr^h 

m = (—-— )m ww 

In each case we can find the functional form of h(r ) that will produce the simplest transport 
and can relate that value to the diffusion constant or the wave velocity. Our example solutions 
include this simple case and a case where the mass and density evolution are more complex. 

The choice hd(r) = ADq/t 5 yields a simple (u,w) diffusion equation for mass. From the 
Ricci relations and the field equations we see that 

m/r 2 = 2 tt h d (p r - p ± ) 

indicating that, with the Ricci relations, mass transport is driven by the pressure anisotropy. 
This also illustrates a drawback of the mass transport ansatz. $n is certainly zero for 
isotropic pressures with a string equation of state [Eci. (IHal) with p± = p r ]. The assumed 
relation, $22 = h d (r)Q n, then says that there is no time variation in the mass for isotropic 
pressures. In general, when the two Ricci components are not related, the mass can vary 
with time, even in the case of isotropic pressures. 

We list the assumptions that lead to matter transport. For diffusion the specific choice 
hd = ADq/t^ is used. 


Assumption 

m = AixD^drp 

$22 = ( 4 :D 0 /r 5 )$ n 

d u $22 = h w (r)§ n 

(d u + a 0 )$22 + h t (r)$u = 0 


Transport 

m — D 0 W^ 2 m = 0, p — D 0 V 2 p = 0 
rh — DQ\7~ 2 m = 0 
fh — (/i„,/4)r 3 V _2 m = 0 
m — (/q/4)rV -2 m + a^iin = 0 
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We have seen that imposing a Ricci relation provides a broad arena to investigate a range 
of atmospheric transport processes in the generalized Vaidya spacetime and is a rich source 


of new analytic mass solutions. The mass solutions that we presented focused on the growth 


or depletion of an atmosphere around a central object. They can be used to describe the 
behavior, for example, of isolated black hole atmospheres but also offer simple models of 
galaxies with a massive black hole at the center. The relationships we imposed on the Ricci 
components were investigated in terms of mass transport although the fundamental rela¬ 
tionships described the evolution of the Ricci components themselves. The Ricci evolution 
is an interesting avenue for further investigation as they are input functions for the Riemann 
invariants. The evolution of the invariants and their syzigies will be discussed elsewhere. 

APPENDIX A: GENERALIZED VAIDYA PRINCIPAL NULL FRAME 

The principal null frame Eq. (TTjl of the Petrov type D metric gf^ obeys 


W = (A' gv /2)l„l„ - {l/r){m^rh v + rh^m v ), 
n^v = ~(A' gv /2)n lJl l„ + {A&/2r)(m li fh v + fh ll m l ,), 
m KV = {A gy /2r)l tl m v - (1 + (cot d/V2r){;m^in v - m^rh u ). 


(Ala) 

(Alb) 

(Ale) 


with both principal null vectors l^ and geodesic. 


For tetrad {u, f, d, (p } and metric gf^ = v^v u — ryfy — dfj,d u — the basis vectors are 

related by 


Vftdx 11 = A l J?du + A^ /2 dr = A^ ,2 [n^ + (A gv /2)l lx ]dx^ : 
f^dx^ = A^ ,2 dr = A~^ 2 [n^ - (A gv /2)/ M ]dx M , 
d^dx^ = rdd = (l/ v /2)(m M + fh^dx^, 

(p^ dx 11 = rsmddcp = —(i/i/2)(m A1 — rh^dx^. 


(A2a) 

(A2b) 

(A2c) 

(A2d) 


The kinematics of the v flow, acceleration, expansion, and shear, are described by 




(A3) 
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where 


af = [fn/r + A gv d r (m / r^A^^f ^, 
cr^ = (0/3)(-2r''f„ + + ^(/y), 

0 = (m/r)A^ /2 . 


(A4a) 

(A4b) 

(A4c) 


Spherical symmetry allows the function m(u,r ) to be identified as the mass within two- 
surfaces of constant u and r, and invariantly defined from the sectional curvature of those 
surfaces: 


2 m/ r 1 2 3 4 5 6 7 8 9 10 11 = Rapp'd 01 (p 13 & tp”. 
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